Cent. Eur. J. Phys. • 1-4 
Author version 



Central European Journal of Physics 



Centrality dependence of freeze-out parameters from 
Au+Au collisions at ^/sjj^j =7.7, 11.5 and 39 GeV 

Research Article 



Lokesh Kumar* (for the STAR collaboration) 

Kent State University, Kent, OH-44242, USA 



Abstract: The RHIC beam energy scan program in its first phase collected data for Au+Au collisions at beam 

energies of 7.7, 11.5 and 39 GeV. The event statistics collected at these lower energies allow us to study 
the centrality dependence of various observables in detail, and compare to fixed-target experiments at 
SPS for similar beam energies. The chemical and kinetic freeze-out parameters can be extracted from the 
experimentally measured yields of identified hadrons within the framework of thermodynamical models. 
These then provide information about the system at the stages of the expansion where inelastic and elastic 
collisions of the constituents cease. We present the centrality dependence of freeze-out parameters for 
Au+Au collisions at midrapidity for -^/Sjvjv = 7.7, 11.5, and 39 GeV from the STAR experiment. The 
chemical freeze-out conditions are obtained by comparing the measured particle ratios (involving n, K, 
p, and p) to those from the statistical thermal model calculations. The kinetic freeze-out conditions are 
extracted at these energies by simultaneously fitting the invariant yields of identified hadrons (n, K, and 
p) using Blast Wave model calculations. 
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1. Introduction 

The experiments at the Relativistic Heavy-Ion Collider (RHIC) aim to look for the signatures of the production 
of a Quark-Gluon Plasma (QGP). The RHIC Beam Energy Scan (BES) [1] program is devoted to study the QCD 
phase diagram which involves searching for the QCD phase boundary and the QCD critical point [2]. The QCD 
phase diagram is usually represented by plotting temperature T vs. baryon chemical potential /j,b- There are two 
main phases expected in the QCD phase diagram - QGP and hadronic gas phase. The plan of the RHIC BES 
is to collect data at different center-of-mass energies and look for the possible signatures of a QGP and QCD 
critical point. It might then also be possible to locate the beam energy below which there is no QGP formation. 
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To access the phase diagram, one needs T and fj,B corresponding to the beam energy, which can be obtained from 
the spectra and ratios of produced particles. 

The bulk properties of particle production in heavy-ion collisions can be studied using identified particle spectra. 
The measurements of particle abundances and transverse momentum distributions could provide information 
about the final stages of the collision evolution at chemical and kinetic freeze-out. The statistical thermal model [3] 
has successfully described the ratios of hadron yields produced in heavy-ion collisions from lower to higher energies. 
The centrality dependence of model parameters, especially at low energies could provide further understanding 
of the collision dynamics. 

In this paper, we present the first measurements of centrality dependence of freeze-out parameters at ^snn =7.7, 
11.5 and 39 GeV. The results are presented for midrapidity \y\ < 0.1 region. The data are taken by the STAR 
experiment for Au+Au collisions in the year 2010. The total events analyzed are about 4 M, 8 M, and 10 M, 
respectively for ^/snn — 7.7, 11.5, and 39 GeV. The centrality selection is done using the uncorrected charged 
track multiplicity measured event-wise in the TPC within < 0.5. Centrality classes represent the fractions of 
this multiplicity distribution. The average numbers of participating nucleons (A r pa rt) are obtained by comparing 
the multiplicity distribution with that from Glauber Monte simulation [1]. 

2. Results and Discussions 

The transverse momentum spectra of hadrons (tt, K, and p) are obtained for the Au+Au collisions at ^snn =7.7, 
11.5 and 39 GeV and different centralities [4]. The proton spectra are not feed-down corrected for weak decays. 
The particle yields and ratios are used to obtain the freeze-out conditions. The statistical thermal model (THER- 
MUS) [3] is used to extract the chemical freeze-out conditions at these energies. Kinetic freeze-out conditions are 
obtained by fitting the spectra with the Blast- Wave (BW) model calculations [5]. 

2.1. Chemical freeze-out 

The scenario after the collisions when the inelastic collisions among the particles stop is called the chemical 
freeze-out. At this stage, the particle yields and ratios are finalized and do not change afterwards. The system 
can be described by the thermal equilibrium model (THERMUS). The model assumes thermal and chemical 
equilibrium. The main fit parameters are chemical freeze-out parameter T c h, baryonic chemical potential fiB, 
strangeness chemical potential [is, and strangeness suppression factor 7 S . The grand-canonical (GC) approach is 
used to fit the experimental particle ratios /ir + , K~ / K + , p/p, K~ /ir~ , and p/n~). For GC ensemble, the 
baryon number, charge number, and strangeness content are conserved on average. The model explains the data 
very well. 

Figure 1 shows the centrality dependence of chemical freeze-out temperature (left panel) and baryonic chemical 
potential (right panel). The T c h increases with increasing energy. It also shows a slight increase as we go from 
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Figure 1. Left: T c h as function of (N paT t) for ^ sn n —7.7, 11.5 and 39 GeV. Right: as function of {A^p ar t) for *J sjy jv —7.7, 
11.5 and 39 GeV. Errors represent the statistical and systematic errors added in quadrature. 



peripheral to central collisions for all energies. The baryonic chemical potential increases with decreasing energy. 
This is because of large baryon stopping at mid-rapidity at low energies. The ^ib also shows a slight increase from 
peripheral to central collisions for these energies. 

Figure 2 (left panel) shows the T c h vs. (j,b plot. Results are shown for Au+Au collisions at ^snn =7.7, 11.5, 39 
GeV and 200 GeV for different centralities. We see that the 200 GeV data do not show any centrality dependence 
for T c h and fiB- The 39 GeV data show a slight dependence of freeze-out parameters on centrality. However, when 
we go towards lower energies (11.5 and 7.7 GeV), we observe a clear centrality dependence of the T c h and /j,b- 
This is the first observation of this behavior by any experiment. It may be noted that the freeze-out parameters 
shown here are obtained only from the ratios involving tt, K,p and p. It will be interesting to see the behavior of 
these parameters after including more ratios involving strange particles such as A and K° in THERMUS fits. 

2.2. Kinetic freeze-out 

Kinetic freeze-out represents the scenario when the interactions among the particles cease. After this time, 
the spectral shapes of the particle species do not change. The kinetic freeze-out conditions are obtained by 
fitting simultaneously the n,K,p spectra with the Blast- Wave model calculations. The BW model describes the 
spectral shapes assuming a locally thermalized source with a common transverse flow velocity field. The main 
fit parameters are the kinetic freeze-out parameter T^ n , the average flow velocity (/3), and the velocity profile 
n. The right panel of Fig. 2 shows the variation of kinetic freeze-out temperature as a function of average flow 
velocity. Results are shown for different energies (7.7, 11.5, 39, 62.4 and 200 GeV) and for different centralities. 
2ki n decreases with increase in energy. It also decreases as we go from peripheral to central collisions. The {j3) 
increases with increase in energy as well as collision centrality. The figure suggests that the higher value of Tki n 
corresponds to the lower value of (/?) and vice-versa. 
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Fi gure 2. Left: Variation of T c h with /ab for different centralities at ^/siviv" —7.7, 11.5, 39 and 200 GeV. Right: Variation of Tki n 
with {/3) for different centralities at V /JKTKT =7.7, 11.5, 39, 62.4 and 200 GeV. The 62.4 and 200 GeV results are taken 
from the Rcf. [6]. Errors represent the statistical and systematic errors added in quadrature. 



3. Summary 

In summary, the centrality dependence of freeze-out parameters (kinetic and chemical) are presented for the BES 
energies 7.7, 11.5 and 39 GeV. The chemical freeze-out conditions are obtained using the particle ratios. For the 
first time, a clear centrality dependence of freeze-out parameters is observed at low energies. New measurements 
have extended the fiB range covered by the RHIC data from 20-400 MeV in the QCD phase diagram. The kinetic 
freeze-out conditions are obtained by using the particle px spectra. The kinetic freeze-out temperature decreases 
with increase in energy and centrality. Average flow velocity increases with increase in energy and centrality. 
We acknowledge the support from DOE and DAE-BRNS project sanction No. 2010/21/15-BRNS/2026. 
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